Elucidating the genetic architecture of DNA methylation to identify promising molecular mechanisms of disease.

DNA methylation commonly occurs at cytosine-phosphate-guanine sites (CpGs) that can serve as biomarkers for many diseases. We analyzed whole genome sequencing data to identify DNA methylation quantitative trait loci (mQTLs) in 4126 Framingham Heart Study participants. Our mQTL mapping identified 94,362,817 cis-mQTLvariant-CpG pairs (for 210,156 unique autosomal CpGs) at P < 1e-7 and 33,572,145 trans-mQTL variant-CpG pairs (for 213,606 unique autosomal CpGs) at P < 1e-14. Using cis-mQTL variants for 1258 CpGs associated with seven cardiovascular disease (CVD) risk factors, we found 104 unique CpGs that colocalized with at least one CVD trait. For example, cg11554650 (PPP1R18) colocalized with type 2 diabetes, and was driven by a single nucleotide polymorphism (rs2516396). We performed Mendelian randomization (MR) analysis and demonstrated 58 putatively causal relations of CVD risk factor-associated CpGs to one or more risk factors (e.g., cg05337441 [APOB] with LDL; MR P = 1.2e-99, and 17 causal associations with coronary artery disease (e.g. cg08129017 [SREBF1] with coronary artery disease; MR P = 5e-13). We also showed that three CpGs, e.g., cg14893161 (PM20D1), are putatively causally associated with COVID-19 severity. To assist in future analyses of the role of DNA methylation in disease pathogenesis, we have posted a comprehensive summary data set in the National Heart, Lung, and Blood Institute's BioData Catalyst.

Age and sex effects on DNA methylation sites linked to genes implicated in severe COVID-19 and SARS-CoV-2 host cell entry.

Male sex and advanced age are associated with severe symptoms of COVID-19. Sex and age also exhibit substantial associations with genome-wide DNA methylation (DNAm) differences in humans. Using a random sample of Illumina EPIC-based genome-wide methylomes from peripheral whole blood of 1,976 parents, participating in The Norwegian Mother, Father and Child Cohort Study (MoBa), we explored whether DNAm in genes linked to SARS-CoV-2 host cell entry and to severe COVID-19 were associated with sex and age. This was carried out by testing 1,572 DNAm sites (CpGs) located near 45 genes for associations with age and sex. We found that DNAm in 281 and 231 of 1,572 CpGs were associated (pFDR<0.01) with sex and aging, respectively. CpGs linked to SARS-CoV-2 host cell entry genes were all associated with age and sex, except for the ACE2 receptor gene (located on the X-chromosome), which was only associated with sex (pFDR<0.01). Furthermore, we examined whether 1,487 autosomal CpGs associated with host-cell entry and severe COVID-19 were more or less associated with sex and age than what would be expected from the same number of randomly sampled genome-wide CpGs. We found that the CpGs associated with host-cell entry and severe COVID-19 were not more or less associated with sex (R2 = 0.77, p = 0.09) than the CpGs sampled from random genomic regions; age was actually found to be significantly less so (R2 = 0.36, p = 0.04). Hence, while we found wide-spread associations between sex and age at CpGs linked to genes implicated with SARS-CoV-2 host cell entry and severe COVID-19, the effect from the sum of these CpGs was not stronger than that from randomly sampled CpGs; for age it was significantly less so. These findings could suggest that advanced age and male sex may not be unsurmountable barriers for the SARS-CoV-2 virus to evolve increased infectiousness.

Persistent blood DNA methylation changes one year after SARS-CoV-2 infection.

We recently reported the COVID-19-induced circulating leukocytes DNA methylation profile. Here, we hypothesized that some of these genes would persist differentially methylated after disease resolution. Fifteen participants previously hospitalized for SARS-CoV-2 infection were epityped one year after discharge. Of the 1505 acute illness-induced differentially methylated regions (DMRs) previously identified, we found 71 regions with persisted differentially methylated, with an average of 7 serial CpG positions per DMR. Sixty-four DMRs persisted hypermethylated, and 7 DMR persisted hypomethylated. These data are the first reported evidence that DNA methylation changes in circulating leukocytes endure long after recovery from acute illness.

S-farnesylation is essential for antiviral activity of the long ZAP isoform against RNA viruses with diverse replication strategies.

The zinc finger antiviral protein (ZAP) is a broad inhibitor of virus replication. Its best-characterized function is to bind CpG dinucleotides present in viral RNAs and, through the recruitment of TRIM25, KHNYN and other cofactors, target them for degradation or prevent their translation. The long and short isoforms of ZAP (ZAP-L and ZAP-S) have different intracellular localization and it is unclear how this regulates their antiviral activity against viruses with different sites of replication. Using ZAP-sensitive and ZAP-insensitive human immunodeficiency virus type I (HIV-1), which transcribe the viral RNA in the nucleus and assemble virions at the plasma membrane, we show that the catalytically inactive poly-ADP-ribose polymerase (PARP) domain in ZAP-L is essential for CpG-specific viral restriction. Mutation of a crucial cysteine in the C-terminal CaaX box that mediates S-farnesylation and, to a lesser extent, the residues in place of the catalytic site triad within the PARP domain, disrupted the activity of ZAP-L. Addition of the CaaX box to ZAP-S partly restored antiviral activity, explaining why ZAP-S lacks antiviral activity for CpG-enriched HIV-1 despite conservation of the RNA-binding domain. Confocal microscopy confirmed the CaaX motif mediated localization of ZAP-L to vesicular structures and enhanced physical association with intracellular membranes. Importantly, the PARP domain and CaaX box together jointly modulate the interaction between ZAP-L and its cofactors TRIM25 and KHNYN, implying that its proper subcellular localisation is required to establish an antiviral complex. The essential contribution of the PARP domain and CaaX box to ZAP-L antiviral activity was further confirmed by inhibition of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication, which replicates in double-membrane vesicles derived from the endoplasmic reticulum. Thus, compartmentalization of ZAP-L on intracellular membranes provides an essential effector function in ZAP-L-mediated antiviral activity against divergent viruses with different subcellular replication sites.

Hydrogel-Based Slow Release of a Receptor-Binding Domain Subunit Vaccine Elicits Neutralizing Antibody Responses Against SARS-CoV-2.

The development of effective vaccines that can be rapidly manufactured and distributed worldwide is necessary to mitigate the devastating health and economic impacts of pandemics like COVID-19. The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein, which mediates host cell entry of the virus, is an appealing antigen for subunit vaccines because it is efficient to manufacture, highly stable, and a target for neutralizing antibodies. Unfortunately, RBD is poorly immunogenic. While most subunit vaccines are commonly formulated with adjuvants to enhance their immunogenicity, clinically-relevant adjuvants Alum, AddaVax, and CpG/Alum are found unable to elicit neutralizing responses following a prime-boost immunization. Here, it has been shown that sustained delivery of an RBD subunit vaccine comprising CpG/Alum adjuvant in an injectable polymer-nanoparticle (PNP) hydrogel elicited potent anti-RBD and anti-spike antibody titers, providing broader protection against SARS-CoV-2 variants of concern compared to bolus administration of the same vaccine and vaccines comprising other clinically-relevant adjuvant systems. Notably, a SARS-CoV-2 spike-pseudotyped lentivirus neutralization assay revealed that hydrogel-based vaccines elicited potent neutralizing responses when bolus vaccines did not. Together, these results suggest that slow delivery of RBD subunit vaccines with PNP hydrogels can significantly enhance the immunogenicity of RBD and induce neutralizing humoral immunity.

PPAR-α Hypermethylation in the Hippocampus of Mice Exposed to Social Isolation Stress Is Associated with Enhanced Neuroinflammation and Aggressive Behavior.

Social behavioral changes, including social isolation or loneliness, increase the risk for stress-related disorders, such as major depressive disorder, posttraumatic stress disorder (PTSD), and suicide, which share a strong neuroinflammatory etiopathogenetic component. The peroxisome-proliferator activated receptor (PPAR)-α, a newly discovered target involved in emotional behavior regulation, is a ligand-activated nuclear receptor and a transcription factor that, following stimulation by endogenous or synthetic ligands, may induce neuroprotective effects by modulating neuroinflammation, and improve anxiety and depression-like behaviors by enhancing neurosteroid biosynthesis. How stress affects epigenetic mechanisms with downstream effects on inflammation and emotional behavior remains poorly understood. We studied the effects of 4-week social isolation, using a mouse model of PTSD/suicide-like behavior, on hippocampal PPAR-α epigenetic modification. Decreased PPAR-α expression in the hippocampus of socially isolated mice was associated with increased levels of methylated cytosines of PPAR-α gene CpG-rich fragments and deficient neurosteroid biosynthesis. This effect was associated with increased histone deacetylases (HDAC)1, methyl-cytosine binding protein (MeCP)2 and decreased ten-eleven translocator (TET)2 expression, which favor hypermethylation. These alterations were associated with increased TLR-4 and pro-inflammatory markers (e.g., TNF-α,), mediated by NF-κB signaling in the hippocampus of aggressive mice. This study contributes the first evidence of stress-induced brain PPAR-α epigenetic regulation. Social isolation stress may constitute a risk factor for inflammatory-based psychiatric disorders associated with neurosteroid deficits, and targeting epigenetic marks linked to PPAR-α downregulation may offer a valid therapeutic approach.

Epigenetic Clocks Are Not Accelerated in COVID-19 Patients.

Age is a major risk factor for severe outcome of the 2019 coronavirus disease (COVID-19). In this study, we followed the hypothesis that particularly patients with accelerated epigenetic age are affected by severe outcomes of COVID-19. We investigated various DNA methylation datasets of blood samples with epigenetic aging signatures and performed targeted bisulfite amplicon sequencing. Overall, epigenetic clocks closely correlated with the chronological age of patients, either with or without acute respiratory distress syndrome. Furthermore, lymphocytes did not reveal significantly accelerated telomere attrition. Thus, these biomarkers cannot reliably predict higher risk for severe COVID-19 infection in elderly patients.

SARS-CoV-2 hot-spot mutations are significantly enriched within inverted repeats and CpG island loci.

SARS-CoV-2 is an intensively investigated virus from the order Nidovirales (Coronaviridae family) that causes COVID-19 disease in humans. Through enormous scientific effort, thousands of viral strains have been sequenced to date, thereby creating a strong background for deep bioinformatics studies of the SARS-CoV-2 genome. In this study, we inspected high-frequency mutations of SARS-CoV-2 and carried out systematic analyses of their overlay with inverted repeat (IR) loci and CpG islands. The main conclusion of our study is that SARS-CoV-2 hot-spot mutations are significantly enriched within both IRs and CpG island loci. This points to their role in genomic instability and may predict further mutational drive of the SARS-CoV-2 genome. Moreover, CpG islands are strongly enriched upstream from viral ORFs and thus could play important roles in transcription and the viral life cycle. We hypothesize that hypermethylation of these loci will decrease the transcription of viral ORFs and could therefore limit the progression of the disease.

Methylation of Host Genes Associated with Coronavirus Infection from Birth to 26 Years.

DNA methylation (DNAm) patterns over time at 1146 CpGs on coronavirus-related genes were assessed to understand whether the varying differences in susceptibility, symptoms, and the outcomes of the SARS-CoV-2 infection in children and young adults could be explained through epigenetic alterations in a host cell's transcriptional apparatus to coronaviruses. DNAm data from the Isle of Wight birth cohort (IOWBC) at birth, 10, 18, and 26 years of age were included. Linear mixed models with repeated measurements stratified by sex were used to examine temporal patterns, and cluster analysis was performed to identify CpGs following similar patterns. CpGs on autosomes and sex chromosomes were analyzed separately. The association of identified CpGs and expression of their genes were evaluated. Pathway enrichment analyses of the genes was conducted at FDR = 0.05. DNAm at 635 of the 1146 CpGs on autosomes showed statistically significant time effects (FDR = 0.05). The 635 CpGs were classified into five clusters with each representing a unique temporal pattern of DNAm. Of the 29 CpGs on sex chromosomes, DNAm at seven CpGs in males and eight CpGs in females showed time effects (FDR = 0.05). Sex-specific and non-specific associations of DNAm with gene expression were found at 24 and 93 CpGs, respectively. Genes which mapped the 643 CpGs represent 460 biological processes. We suggest that the observed variability in DNAm with advancing age may partially explain differing susceptibility, disease severity, and mortality of coronavirus infections among different age groups.

An epigenome-wide DNA methylation study of patients with COVID-19.

In the early 2000s, emerging SARS-CoV-2, which is highly pathogenic, posed a great threat to public health. During COVID-19, epigenetic regulation is deemed to be an important part of the pathophysiology and illness severity. Using the Illumina Infinium Methylation EPIC BeadChip (850 K), we investigated genome-wide differences in DNA methylation between healthy subjects and COVID-19 patients with different disease severities. We conducted a combined analysis and selected 35 "marker" genes that could indicate a SARS-CoV-2 infection, including 12 (ATHL1, CHN2, CHST15, CPLX2, CRHR2, DCAKD, GNAI2, HECW1, HYAL1, MIR510, PDE11A, and SMG6) situated in the promoter region. The functions and pathways of differentially methylated genes were enriched in biological processes, signal transduction, and the immune system. In the "Severe versus Mild" group, differentially methylated genes, after eliminating duplicates, were used for PPI analyses. The four hub genes (GNG7, GNAS, PRKCZ, and PRKAG2) that had the highest degree of nodes were identified and among them, GNG7 and GNAS genes expressions were also downregulated in the severe group in sequencing results. Above all, the results suggest that GNG7 and GNAS may play a non-ignorable role in the progression of COVID-19. In conclusion, the identified key genes and related pathways in the current study can be used to study the molecular mechanisms of COVID-19 and may provide possibilities for specific treatments.

Evidence for Biological Age Acceleration and Telomere Shortening in COVID-19 Survivors.

The SARS-CoV-2 infection determines the COVID-19 syndrome characterized, in the worst cases, by severe respiratory distress, pulmonary and cardiac fibrosis, inflammatory cytokine release, and immunosuppression. This condition has led to the death of about 2.15% of the total infected world population so far. Among survivors, the presence of the so-called persistent post-COVID-19 syndrome (PPCS) is a common finding. In COVID-19 survivors, PPCS presents one or more symptoms: fatigue, dyspnea, memory loss, sleep disorders, and difficulty concentrating. In this study, a cohort of 117 COVID-19 survivors (post-COVID-19) and 144 non-infected volunteers (COVID-19-free) was analyzed using pyrosequencing of defined CpG islands previously identified as suitable for biological age determination. The results show a consistent biological age increase in the post-COVID-19 population, determining a DeltaAge acceleration of 10.45 ± 7.29 years (+5.25 years above the range of normality) compared with 3.68 ± 8.17 years for the COVID-19-free population (p < 0.0001). A significant telomere shortening parallels this finding in the post-COVID-19 cohort compared with COVID-19-free subjects (p < 0.0001). Additionally, ACE2 expression was decreased in post-COVID-19 patients, compared with the COVID-19-free population, while DPP-4 did not change. In light of these observations, we hypothesize that some epigenetic alterations are associated with the post-COVID-19 condition, particularly in younger patients (< 60 years).

SARS-CoV-2 and other human coronavirus show genome patterns previously associated to reduced viral recognition and altered immune response.

A new pandemic caused by the betacoronavirus SARS-CoV-2 originated in China in late 2019. Although often asymptomatic, a relevant percentage of affected people can develop severe pneumonia. Initial evidence suggests that dysregulation of the immune response could contribute to the pathogenesis, as previously demonstrated for SARS-CoV. The presence of genome composition features involved in delaying viral recognition is herein investigated for human coronaviruses (HCoVs), with a special emphasis on SARS-CoV-2. A broad collection of HCoVs polyprotein, envelope, matrix, nucleocapsid and spike coding sequences was downloaded and several statistics representative of genome composition and codon bias were investigated. A model able to evaluate and test the presence of a significant under- or over-representation of dinucleotide pairs while accounting for the underlying codon bias and protein sequence was also implemented. The study revealed the significant under-representation of CpG dinucleotide pair in all HcoV, but especially in SARS-CoV and even more in SARS-CoV-2. The presence of forces acting to minimize CpG content was confirmed by relative synonymous codon usage pattern. Codons containing the CpG pair were severely under-represented, primarily in the polyprotein and spike coding sequences of SARS-CoV-2. Additionally, a significant under-representation of the TpA pair was observed in the N and S region of SARS-CoV and SARS-CoV-2. Increasing experimental evidence has proven that CpG and TpA are targeted by innate antiviral host defences, contributing both to RNA degradation and RIG-1 mediated interferon production. The low content of these dinucleotides could contribute to a delayed interferon production, dysregulated immune response, higher viral replication and poor outcome. Significantly, the RIG-1 signalling pathway was proven to be defective in elderlies, suggesting a likely interaction between limited viral recognition and lower responsiveness in interferon production that could justify the higher disease severity and mortality in older patients.

Epigenome-wide association study of COVID-19 severity with respiratory failure.

BACKGROUND: Patients infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the coronavirus disease 2019 (COVID-19), exhibit a wide spectrum of disease behaviour. Since DNA methylation has been implicated in the regulation of viral infections and the immune system, we performed an epigenome-wide association study (EWAS) to identify candidate loci regulated by this epigenetic mark that could be involved in the onset of COVID-19 in patients without comorbidities. METHODS: Peripheral blood samples were obtained from 407 confirmed COVID-19 patients ≤ 61 years of age and without comorbidities, 194 (47.7%) of whom had mild symptomatology that did not involve hospitalization and 213 (52.3%) had a severe clinical course that required respiratory support. The set of cases was divided into discovery (n = 207) and validation (n = 200) cohorts, balanced for age and sex of individuals. We analysed the DNA methylation status of 850,000 CpG sites in these patients. FINDINGS: The DNA methylation status of 44 CpG sites was associated with the clinical severity of COVID-19. Of these loci, 23 (52.3%) were located in 20 annotated coding genes. These genes, such as the inflammasome component Absent in Melanoma 2 (AIM2) and the Major Histocompatibility Complex, class I C (HLA-C) candidates, were mainly involved in the response of interferon to viral infection. We used the EWAS-identified sites to establish a DNA methylation signature (EPICOVID) that is associated with the severity of the disease. INTERPRETATION: We identified DNA methylation sites as epigenetic susceptibility loci for respiratory failure in COVID-19 patients. These candidate biomarkers, combined with other clinical, cellular and genetic factors, could be useful in the clinical stratification and management of patients infected with the SARS-CoV-2. FUNDING: The Unstoppable campaign of the Josep Carreras Leukaemia Foundation, the Cellex Foundation and the CERCA Programme/Generalitat de Catalunya.

The Heterogeneous Landscape and Early Evolution of Pathogen-Associated CpG Dinucleotides in SARS-CoV-2.

COVID-19 can lead to acute respiratory syndrome, which can be due to dysregulated immune signaling. We analyze the distribution of CpG dinucleotides, a pathogen-associated molecular pattern, in the SARS-CoV-2 genome. We characterize CpG content by a CpG force that accounts for statistical constraints acting on the genome at the nucleotidic and amino acid levels. The CpG force, as the CpG content, is overall low compared with other pathogenic betacoronaviruses; however, it widely fluctuates along the genome, with a particularly low value, comparable with the circulating seasonal HKU1, in the spike coding region and a greater value, comparable with SARS and MERS, in the highly expressed nucleocapside coding region (N ORF), whose transcripts are relatively abundant in the cytoplasm of infected cells and present in the 3'UTRs of all subgenomic RNA. This dual nature of CpG content could confer to SARS-CoV-2 the ability to avoid triggering pattern recognition receptors upon entry, while eliciting a stronger response during replication. We then investigate the evolution of synonymous mutations since the outbreak of the COVID-19 pandemic, finding a signature of CpG loss in regions with a greater CpG force. Sequence motifs preceding the CpG-loss-associated loci in the N ORF match recently identified binding patterns of the zinc finger antiviral protein. Using a model of the viral gene evolution under human host pressure, we find that synonymous mutations seem driven in the SARS-CoV-2 genome, and particularly in the N ORF, by the viral codon bias, the transition-transversion bias, and the pressure to lower CpG content.

Evolution of Transmissible Gastroenteritis Virus (TGEV): A Codon Usage Perspective.

Transmissible gastroenteritis virus (TGEV) is a coronavirus associated with diarrhea and high mortality in piglets. To gain insight into the evolution and adaptation of TGEV, a comprehensive analysis of phylogeny and codon usage bias was performed. The phylogenetic analyses of maximum likelihood and Bayesian inference displayed two distinct genotypes: genotypes I and II, and genotype I was classified into subtypes Ia and Ib. The compositional properties revealed that the coding sequence contained a higher number of A/U nucleotides than G/C nucleotides, and that the synonymous codon third position was A/U-enriched. The principal component analysis based on the values of relative synonymous codon usage (RSCU) showed the genotype-specific codon usage patterns. The effective number of codons (ENC) indicated moderate codon usage bias in the TGEV genome. Dinucleotide analysis showed that CpA and UpG were over-represented and CpG was under-represented in the coding sequence of the TGEV genome. The analyses of Parity Rule 2 plot, ENC-plot, and neutrality plot displayed that natural selection was the dominant evolutionary driving force in shaping codon usage preference in genotypes Ia and II. In addition, natural selection played a major role, while mutation pressure had a minor role in driving the codon usage bias in genotype Ib. The codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses suggested that genotype I might be more adaptive to pigs than genotype II. Current findings contribute to understanding the evolution and adaptation of TGEV.

SARS-CoV-2 Is Restricted by Zinc Finger Antiviral Protein despite Preadaptation to the Low-CpG Environment in Humans.

Recent evidence shows that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is sensitive to interferons (IFNs). However, the most effective types of IFNs and the underlying antiviral effectors remain to be defined. Here, we show that zinc finger antiviral protein (ZAP), which preferentially targets CpG dinucleotides in viral RNA sequences, restricts SARS-CoV-2. We further demonstrate that ZAP and its cofactors KHNYN and TRIM25 are expressed in human lung cells. Type I, II, and III IFNs all strongly inhibited SARS-CoV-2 and further induced ZAP expression. Comprehensive sequence analyses revealed that SARS-CoV-2 and its closest relatives from horseshoe bats showed the strongest CpG suppression among all known human and bat coronaviruses, respectively. Nevertheless, endogenous ZAP expression restricted SARS-CoV-2 replication in human lung cells, particularly upon treatment with IFN-α or IFN-γ. Both the long and the short isoforms of human ZAP reduced SARS-CoV-2 RNA expression levels, but the former did so with greater efficiency. Finally, we show that the ability to restrict SARS-CoV-2 is conserved in ZAP orthologues of the reservoir bat and potential intermediate pangolin hosts of human coronaviruses. Altogether, our results show that ZAP is an important effector of the innate response against SARS-CoV-2, although this pandemic pathogen emerged from zoonosis of a coronavirus that was preadapted to the low-CpG environment in humans.IMPORTANCE Although interferons inhibit SARS-CoV-2 and have been evaluated for treatment of coronavirus disease 2019 (COVID-19), the most effective types and antiviral effectors remain to be defined. Here, we show that IFN-γ is particularly potent in restricting SARS-CoV-2 and in inducing expression of the antiviral factor ZAP in human lung cells. Knockdown experiments revealed that endogenous ZAP significantly restricts SARS-CoV-2. We further show that CpG dinucleotides which are specifically targeted by ZAP are strongly suppressed in the SARS-CoV-2 genome and that the two closest horseshoe bat relatives of SARS-CoV-2 show the lowest genomic CpG content of all coronavirus sequences available from this reservoir host. Nonetheless, both the short and long isoforms of human ZAP reduced SARS-CoV-2 RNA levels, and this activity was conserved in horseshoe bat and pangolin ZAP orthologues. Our findings indicating that type II interferon is particularly efficient against SARS-CoV-2 and that ZAP restricts this pandemic viral pathogen might promote the development of effective immune therapies against COVID-19.

Viral replication of SARS-CoV-2 could be self-limitative - The role of the renin-angiotensin system on COVID-19 pathophysiology.

Currently, the world is suffering with one of the biggest pandemics of recent history. Caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the coronavirus disease 2019 (COVID-19) is provoking devastating consequences on economic and social fields throughout all continents. Therefore, pathophysiological knowledge about COVID-19 is imperative for better planning of preventive measures, diagnosis, and therapeutics of the disease. Based on previous studies, this work proposes new hypothesis related to the role of the renin-angiotensin system on the pathophysiology of COVID-19, and its purpose is to enrich the discussion and to offer alternative ways for experimental and clinical studies aiming at the formulation of new diagnosis and/or treatment methods.

Human SARS-CoV-2 has evolved to reduce CG dinucleotide in its open reading frames.

The outbreak of COVID-19 has brought great threat to human health. Its causative agent is a severe acute respiratory syndrome-related coronavirus which has been officially named SARS-CoV-2. Here we report the discovery of extremely low CG abundance in its open reading frames. We found that CG reduction in SARS-CoV-2 is achieved mainly through mutating C/G into A/T, and CG is the best target for mutation. Meanwhile, 5'-untranslated region of SARS-CoV-2 has high CG content and is capable of forming an internal ribosome entry site (IRES) to recruit host ribosome for translating its RNA. These features allow SARS-CoV-2 to reproduce efficiently in host cells, because less energy is consumed in disrupting the stem-loops formed by its genomic RNA. Notably, genomes of cellular organisms also have very low CG abundance, suggesting that mutating C/G into A/T occurs universally in all life forms. Moreover, CG is the dinucleotide related to CpG island, mutational hotspot and single nucleotide polymorphism in cellular organisms. The relationship between these features is worthy of further investigations.

Viral CpG Deficiency Provides No Evidence That Dogs Were Intermediate Hosts for SARS-CoV-2.

Due to the scope and impact of the COVID-19 pandemic there exists a strong desire to understand where the SARS-CoV-2 virus came from and how it jumped species boundaries to humans. Molecular evolutionary analyses can trace viral origins by establishing relatedness and divergence times of viruses and identifying past selective pressures. However, we must uphold rigorous standards of inference and interpretation on this topic because of the ramifications of being wrong. Here, we dispute the conclusions of Xia (2020. Extreme genomic CpG deficiency in SARS-CoV-2 and evasion of host antiviral defense. Mol Biol Evol. doi:10.1093/molbev/masa095) that dogs are a likely intermediate host of a SARS-CoV-2 ancestor. We highlight major flaws in Xia's inference process and his analysis of CpG deficiencies, and conclude that there is no direct evidence for the role of dogs as intermediate hosts. Bats and pangolins currently have the greatest support as ancestral hosts of SARS-CoV-2, with the strong caveat that sampling of wildlife species for coronaviruses has been limited.